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Abstract 

Dynamical models of prototype gravastars made of phantom energy are constructed, in which 
an infinitely thin spherical shell of a perfect fluid with the equation of state p = (1 — 7)0' divides 
the whole spacetime into two regions, the internal region filled with a dark energy (or phantom) 
fluid, and the external Schwarzschild region. It is found that in some cases the models represent 
the "bounded excursion" stable gravastars, where the thin shell is oscillating between two finite 
radii, while in other cases they collapse until the formation of black holes or normal stars. In the 
phase space, the region for the "bounded excursion" gravastars is very small in comparison to that 
of black holes, but not empty, as found in our previous papers. Therefore, although the existence 
of gravastars can not be completely excluded from current analysis, the opposite is not possible 
either, that is, even if gravastars exist, they do not exclude the existence of black holes. 
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I. INTRODUCTION 

As alternatives to black holes, gravastars have received some attention recently par- 
tially due to the tight connection between the cosmological constant and a currently accel- 
erating universe although very strict observational constraints on the existence of such 
stars may exist [3j. ^ 

The pioneer model of gravastar was proposed by Mazur and Mottola (MM) j4] , consisting 
of five layers: an internal core < r < ri , described by the de Sitter universe, an inter- 
mediate thin layer of stiff fluid ri < r < r2 , an external region r > r2 , described by the 
Schwarzschild solution, and two infinitely thin shells, appearing, respectively, on the hyper- 
surfaces r = ri and r = r2. The intermediate layer is constructed in such way that ri is inner 
than the de Sitter horizon, while r2 is outer than the Schwarzschild horizon, eliminating the 
apparent horizon. Configurations with a de Sitter interior have long history which we can 
find, for example, in the work of Dymnikova and Galaktionov j^. After this work, Visser 
and Wiltshire 6|] pointed out that there are two different types of stable gravastars which are 
stable gravastars and "bounded excursion" gravastars. In the spherically symmetric case, 
the motion of the surface of the gravastar can be written in the form 6|, 

la' + Via) = 0, (1) 

where a denotes the radius of the star, and a = da/dr, with r being the proper time of the 
surface. Depending on the properties of the potential V{a), the two kinds of gravastars are 
defined as follows. 

Stable gravastars: In this case, there must exist a radius ao such that 

V{ao) = 0, V'{ao)=0, (ao) > 0, (2) 

where a prime denotes the ordinary differentiation with respect to the indicated argument. 
If and only if there exists such a radius for which the above conditions are satisfied, the 
model is said to be stable. Among other things, VW found that there are many equations 
of state for which the gravastar configurations are stable, while others are not Carter 
studied the same problem and found new equations of state for which the gravastar is stable 



while De Benedictis et al iSIl and Chirenti and RezzoUa 



^ investigated the stability 



of the original model of Mazur and Mottola against axial-perturbations, and found that 
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gravastars are stable to these perturbations, too. Chirenti and RezzoUa also showed that 
their quasi-normal modes differ from those of black holes with the same mass, and thus can 
be used to discern a gravastar from a black hole. 

"Bounded excursion" gravastars: As VW noticed, there is a less stringent notion of 
stability, the so-called "bounded excursion" models, in which there exist two radii Oi and 02 
such that 

F(ai) = 0, F'(ai)<0, l^(a2)=0, l^' (aa) > 0, (3) 

with V{a) < for a G (ai, 02), where 02 > ai. 

Lately, we studied both types of gravastars 0, and found that, such configurations 
can indeed be constructed, although the region for the formation of them is very small in 
comparison to that of black holes. 

Based on the discussions about the gravastar picture some authors have proposed al- 



te.native .nodel. M- A„»ng then,, we can find a Chap,yg.n da* .tar M, a g.ava.ta. 
supported by non-linear electrodynamics [IJ], a gravastar with continuous anisotropic pres- 



sure [l5| and recently, Dzhunushaliev et al. worked on spherically symmetric configurations 
of a phantom scalar field and they found something like a gravastar but ,t was unstableQ. 

In addition, Lobo [17] studied two models for a dark energy fluid. One of them describes 
a homogeneous energy density and the other describes an ad-hoc monotonically decreasing 
energy density, although both of them are with anisotropic pressure. In order to match an 
exterior Schwarzschild spacetime he introduced a thin shell between the interior and the 
exterior spacetimes. 



In this paper, we generalize our previous works [10|, to the case where the equation 
of state of the infinitely thin shell is given by p = {1 — 7)0" with 7 being a constant, the 



interior consists of a phantom energy fluid 17||, while the exterior is still the Schwarzschild 



space. We shall first construct three-layer dynamical models, and then show both types of 
gravastars and black holes exist for various situations. The rest of the paper is organized as 
follows: In Sec. II we present the metrics of the interior and exterior spacetimes, and write 
down the motion of the thin shell in the form of Eq.([T]). In Sec. Ill we show the definitions 
of dark and phantom energy, for the interior fluid and for the shell. In Sec. IV we discuss 
the formation of black holes from standard or phantom energy. In Sec. V we analyze the 
formation of gravastar or normal star from standard or phantom energy. In Sec. VI we 
study special cases where we can not have the "bounded excursion". Finally, in Sec. VII we 



present our conclusions. 



II. DYNAMICAL THREE-LAYER PROTOTYPE GRAVASTARS 



The interior fluid is made of an anisotropic dark energy fluid with a metric given by 17 1 

dsl_ = -fidt^ + f^dr^ + r^d^^, (4) 
where dVL^ = dO"^ + sin^(6')ci(/)^ and 



2\u) 



= (1 + 6r2) — (1 + 26r 

where is a constant, and its physical meaning can be seen from the following equation (I6l). 
Since the mass is given by m(r) = 6r'^/[2(l + 26r^)] then we have that 6 > 0. The corre- 
sponding energy density p, radial and tangential pressures Pr and pt are given, respectively, 
by 

fujh\ ( 3 + 26r2 
^ _ _ /_^\ / a;(3 + 2&r^) \ ^ feV^ 



vry V (1 + 2fer2)2 y 327r [(1 + 26r2)3(l + 6r2)] 
|(1 + cj)(3 + 26r2) [(1 + 3cj) + 2hr'^{\ + u;)] 

-8Lu(5 + 26r2)(l + 6r2)|. (6) 

The exterior spacetime is given by the Schwarzschild metric 

ds\ = -fdv^ + r^dr^ + r^dn^, (7) 

where / = 1 — 2m/ r. The metric of the hypersurface on the shell is given by 

dsl = -dT'' + R'{T)d^f. (8) 

Since ds']_ = ds"^ = we flnd that rs = = R, and 

/it' - f2R' = 1, (9) 
f^"-f = 1' (10) 



5 



where the dot denotes the ordinary differentiation with respect to the proper time. On the 
other hand, the interior and exterior normal vectors to the thin shell are given by 



< = i-R,v,0,0). 



(11) 



Then, the interior and exterior extrinsic curvature are given by 



:i + 2bR'^rVl + bRWP - (1 + 2bR'^rVl + bRH^ 



{2bR^uj + 2bR^ + 3uj + 1)- 



2(1 + bR^y/'^il + 2bR^)R^^ (1 + 26i?2)-2(i + bR^yHR + Rt - Rt, 

i{l + bR^)R 
l + 26i?2 ' 

v{4m^v^ - AmRv^ - 3R^R^ + RH'^){2m - Ry^mR'^ + Rv - Rv, 
-v{2m - R), 



Since [l8| 

[Kee] = K+ - K, 
where M is the mass of the shell, we find that 



-M, 



M = v{2m -R) + 



i{l + bR^)R 



1 + 26i?2 

Then, substituting equations ([9]) and (fTO!) into (fT9!) we get 



(12) 

(13) 
(14) 
(15) 
(16) 
(17) 

(18) 
(19) 



M = -R{l- — +R^ 
R 



1/2 



+ R- 





1 + 


- + 2bR^) 


1/2 


(1 + 6i?2)-K 


^i)/4(l + 26i?2)( 





(20) 



In order to keep the ideas of MM as much as possible, we consider the thin shell as consisting 
of a fluid with the equation of state, p = {1 — 7)0", where a and p denote, respectively, the 
surface energy density and pressure of the shell and 7 is a constant. Then, the equation of 
motion of the shell is given by |l8| 



M + MRRp = ATTR'[T^pu''n'^] = ttR' [T+ulnl - 



-a/3 



u n_ 



(21) 
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where is the four-velocity. Since the interior fluid is made of an anisotropic fluid and the 
exterior is vacuum, we get 

M + 87rRR{l--f)a = 0. (22) 
Recall that a = M/{A-kR'^), we flnd that Eq.([22]) has the solution 



M = kR^^^-^\ 



(23) 



where k is an integration constant. Substituting Eq. (l23i) into Eg. (1201) . and rescaling m, b 
and R as, 



m — s> mk 27-3^ 
b bk^, 
R Rk~^, 

we flnd that it can be written in the form of Eq. ([1]) with a replaced by R, and 

1 



(24) 



V{R, m, Lu, b, 7) 



2R%. 



6^"+') - b'l 



^('^+2)^:^4(7-1) ^('^+l)/2 



(25) 



2^(3.+4)/2^2{,-l)^(..+l)/4 U-.)^(.+l)/2^2 _ ^^.+1) /2 ^2 



-24"")6S"+')/'mi? + biR^ + b2R^ + 2b2mR + 62^?'^^"'^ 



1/2 



R-" - 2b 



{w+2) 



mRb\ 



iui+l)/2 



+2b^t^'^mR + 42^+3) i?4(,-l) _ ^(..+2)^2 ^ ^(..+l)^(^+3)/2^2| _ ^26) 



where 



61 = 1 + bR^, b2 = l + 2bR^ 



(27) 



Clearly, for any given constants m, uj,b and 7, equation (!26ll uniquely determines the collapse 
of the prototype gravastar. Depending on the initial value i?o, the collapse can form either 
a black hole, a gravastar, a Minkowski, or a spacetime flUed with phantom fluid. In the 
last case, the thin shell flrst collapses to a flnite non-zero minimal radius and then expands 
to inflnity. To guarantee that initially the spacetime does not have any kind of horizons, 
cosmo logical or event, we must restrict Rq to the range. 



Rq > 2m, 



(28) 
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where Rq is the initial collapse radius. When m = = b, the thin shell disappears, and the 
whole spacetime is Minkowski. So, in the following we shall not consider this case. 

Since the potential (12^ is so complicated, it is too difficult to study it analytically. 
Instead, in the following we shall study it numerically. Before doing so, we shall 



III. CLASSIFICATIONS OF MATTER, DARK ENERGY, AND PHANTOM EN- 
ERGY FOR ANISOTROPIC FLUIDS 



Recently 20|], the classification of matter, dark and phantom energy for an anisotropic 
fluid was given in terms of the energy conditions. Such a classification is necessary for systems 
where anisotropy is important, and the pressure components may play very important roles 
and can have quite different contributions. In this paper, we will use this classification to 
study the collapse of the dynamical prototype gravastars, constructed in the last section. 
In particular, we define dark energy as a fluid which violates the strong energy condition 
(SEC). From the Raychaudhuri equation, we can see that such defined dark energy always 
exerts divergent forces on time-like or null geodesies. On the other hand, we define phantom 
energy as a fluid that violates at least one of the null energy conditions (NEC's). We shall 
further distinguish phantom energy that satisfies the SEC from that which does not satisfy 
the SEC. We call the former attractive phantom energy, and the latter repulsive phantom 
energy. Such a classification is summarized in Table I. 

For the sake of completeness, in Table II we apply it to the matter field located on the 
thin shell, while in Table III we combine all the results of Tables I and II, and present all 
the possibilities. 

In order to consider the equations (jlj) and ([6]) for describing dark energy stars we must 
analyze carefully the ranges of the parameter uj that in fact furnish the expected fluids. It 
can be shown that the condition p + > is violated for uj < —1 and fulffiled for > — 1, 
for any values of R and b. The conditions p + pt > and p + Pr + 2pt > are satisfied for 
uj < —1 and — l/3<cij<0, for any values of R and b. For the other intervals of u the energy 
conditions depend on very complicated relations of R and b. See figure [TJ This provides 
an explicit example, in which the definition of dark energy must be dealed with great care. 



Another case was provided in a previous work 



Note that in the paper ITj where the 



solution is used for the ffist time to model dark energy star, the author presented a particular 
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TABLE I: This table summarizes the classification of the interior matter field, based on the energy 
conditions 2], where we assume that p > 0. 



Matter 






Condition 1 




Condit 


ion 2 


Condition 3 


Normal Matter 




p 


+ Pr 


+ 




> 





P + Pr 


> 





P + Pt 


> 





Dark Energy 




p 


+ Pr 


+ 


2pt 


< 





P + Pr 


> 





P + Pt 


> 





Repulsive Phantom Energ 


sy 


p 


+ Pr 


+ 


2pt 


< 





P + Pr 


< 





P + Pt 


> 





Repulsive Phantom Energ 


sy 


p 


+ Pr 


+ 


2pt 


< 





P + Pr 


> 





P + Pt 


< 





Repulsive Phantom Energ 


sy 


p 


+ Pr 


+ 


2pt 


< 





P + Pr 


< 





P + Pt 


< 





Attractive Phantom Enerj 


?y 


p 


+ Pr 


+ 


2pt 


> 





P + Pr 


< 





P + Pt 


> 





Attractive Phantom Enerj 


?y 


p 


+ Pr 


+ 


2pt 


> 





P + Pr 


> 





P + Pt 


< 





Attractive Phantom Enerj 


?y 


p 


+ Pr 


+ 


2pt 


> 





P + Pr 


< 





P + Pt 


< 






TABLE II: This table summarizes the classification of matter on the thin shell, based on the 
we assume that p >0. 



Matter 



Condition 1 



Condition 2 



7 



Normal Matter 
Dark Energy 
Repulsive Phantom Energy 
Attractive Phantom Energy 



o- + 2p > 
a + 2p<0 
a + 2p<0 
a + 2p>0 



a + p>0 
a + p>0 
a +p <0 
a +p <0 



-1 or 
7/4 
3 

Not possible 



case for uj = —0.5. However, it is easy to see graphically that for this value of oo the weak 
and strong energy conditions are satisfied for any positive value of the parameter b and the 
coordinate R. Thus, the corresponding solution does not represent a dark energy star as it 
is claimed by the author. Taking several values of u in the intervals — 1 < u; < — 1/3 and 
u; > 0, we could not found any case where the interior dark energy exist. 

In order to fulfill the energy condition a + 2p > of the shell and assuming that p = 
(1 — 7)cr we must have 7 < 3/2. On the other hand, in order to satisfy the condition 
cr + p > 0, we obtain 7 < 2. Hereinafter, we will use only some particular values of the 
parameter 7 which are analyzed in this work. See Table II. 
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P+Pr < 




P+Pr > 




P+I\ >o 


P+Pr > 




P+Pr > 


P+Pr+2Pt > 


P+Pr+2Pt >0 


attractive phantom energy 




standard 






energy 





-1 -1/3 



FIG. 1: In this figure we show the intervals of u for which the weak and strong energy conditions 
are independent of the coordinate R and the parameter b. The condition p + Pr > is violated 
for w < — 1 and fulfilled for w > —1, for any values of R and b. The conditions p + pt > and 

P + Pr + '^Pt > arc satisfied for cu < —1 and — 1/3 < a; < 0, for any values of R and b. For the 
others intervals of uj the analysis of the energy conditions depends on a complex relation of R and 
b. 

In the next sections we will discuss three physical possibilities for the type of system that 
can be formed from the study of the potential V{R, m, lu, b, 7): (a) Black hole or dispersion 
of the matter, (b) Gravastar or normal star and (c) Black hole or phantom gravastar. 

IV. BLACK HOLE OR DISPERSION OF THE MATTER 

For m > rric the potential V{R) is strictly negative as shown in figures 4, 5, 8, 9, 12, 
13, 16, 17, 20, 21, 24, 25, 28, 29, 32, 33, 36, 37, 38 and 39. Then, the collapse always 
forms black holes. For m = ttIc, there are two different possibilities, depending on the choice 
of the initial radius Rq. In particular, if the star begins to collapse with Rq > Rc, it will 
approach to the minimal radius Rc- Once it reaches this point, the shell will stop collapsing. 
However, this point is unstable and any small perturbations will lead the star either to 
expand for ever and leave behind a flat spacetime, or to collapse until R = 0, whereby a 
Schwarzschild black hole is finally formed. On the other hand, if the star begins to collapse 
with 2mc < Rq < Rc as shown in these figures, the star will collapse until a black hole is 
formed. For m < rric, the potentials V{R) for each case have a positive maximal, and the 
equation V{R, m < nic) — has two positive roots Ri^2 with R2 > Ri > 0. There are also 
two possibilities here, depending on the choice of the initial radius Rq. If Rq > R2, the 
star will first contract to its minimal radius R — R2 and then expand to infinity, whereby a 
Minkowski spacetime is finally formed. If 2m < Rq < Ri, the star will collapse continuously 
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until i? = 0, and a black hole will be finally formed. 



V. GRAVASTAR OR NORMAL STAR 



In this case the potential takes the shape given by figures 2, 3, 6, 7, 10, 11, 14, 15, 18, 19, 
22, 23, 26, 27, 30, 31, 34 and 35, from which we can see that V{R) = now can have one, two 
or three real roots, depending on the mass of the shell. For m > rric we have, say, where 
> Ri. If we choose Rq > R3 (for m = rric we have R2 = R3), then the star will not be 
allowed in this region because the potential is greater than the zero. However, if we choose 
Ri < Rq < R2, the collapse will bounce back and forth between R = Ri and R = R2- Such a 
possibility is shown in these figures. This is exactly the so-called "bounded excursion" model 
mentioned in joj], and studied in some details in 10, 11|. Of course, in a realistic situation, the 
star will emit both gravitational waves and particles, and the potential shall be self-adjusted 



to produce a minimum aX R = Rstatic where V (R = R 



static J 



= V' {R = Rstatic) whereby 

a gravastar or a normal star is finally formed [6|, 110|, lUl , although in 10, lUj the potential 
tends to —00 when R tends to 00. Here it is completely different since the potential now 
tends to +00 when R tends to 00. Thus, in the cases studied here we do not have situations 



where the star expands leaving behind a fiat spacetime, as in 10|, lll |. 



VI. BLACK HOLE OR PHANTOM GRAVASTAR 



In this case the potential takes the shape given by figures 40, 41, 42 and 43, from which 
we can see that V{R) = now has four real roots, say, Ri, where > -R,. If we choose 
Rq > R4,, then again the star will not be allowed in this region because the potential is 
greater than zero. However, if we choose R3 < Rq < R4, the collapse will bounce back and 
forth between R = R^ and R = R^, as in the previous case. But, if we choose R2 < Rq < R3, 
we can note that this region is forbidden because either the potential is imaginary or greater 
than zero. However, if we choose Ri < Rq < R2, the collapse will bounce back and forth 
between R = Ri and R = R2- U Rq < Ri the system will collapse until R = 0, whereby a 
Schwarzschild black hole is finally formed. 
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TABLE III: This table summarizes all possible kind of energy of the interior fluid and of the shell. 
The boldface figure numbers represent stable structures. 



3ase 


Interior Energy 


Shell Energy 


Figures 


Structures 


A 


Standard 


Standard 


2,3,6,7,10, 11,14,15,18,19 


Normal Star 


B 


Standard 


Dark 


4,8,12,16,20 


Black Hole/Dispersion 


C 


Standard 


Repulsive Phantom 


5,9,13,17,21 


Black Hole/Dispersion 


D 


Dark 


Standard 




Interior not found 


E 


Dark 


Dark 




Interior not found 


F 


Dark 


Repulsive Phantom 




Interior not found 


G 


Repulsive Phantom 


Standard 


22,23,26,27, 30,31,34,35 


Gravastar 


H 


Repulsive Phantom 


Dark 


24,28,32,36 


Black Hole/Dispersion 


I 


Repulsive Phantom Repulsive Phantom 


25,29,33,37 


Black Hole/Dispersion 


J 


Attractive Phantom 


Standard 


41, 43, 40, 42 


Gravastar or Black Hole 


K 


Attractive Phantom 


Dark 


38 


Black Hole/Dispersion 


L 


Attractive Phantom 


Repulsive Phantom 


39 


Black Hole/Dispersion 



VII. CONCLUSIONS 



In this paper, we have studied the problem of the stability of gravastars by constructing 
dynamical three- layer models of VW 6^, which consists of an internal phantom fluid, a 
dynamical infinitely thin shell of perfect fluid with the equation of state p = {1 — 7)0", and 
an external Schwarzschild space. 

It must be noted that, although phantoms have been used widely to explain the late 
cosmic acceleration of the universe, microscopic models of them, either for an anisotropic or 
a perfect fluid, has not been constructed, yet, if there is any. 

We have shown explicitly that the final output can be a black hole, a "bounded excursion" 
stable gravastar, a Minkowski, or a phantom spacetime, depending on the total mass m of 
the system, the parameter u, the constant b, the parameter 7 and the initial position Ro of 
the dynamical shell. All these possibilities have non-zero measurements in the phase space 
of m, b, uj, 7 and Rq, although the region of gravastars is very small in comparison with that 
of black holes. All the results can be summarized in Table III. An interesting result that we 
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FIG. 2: The potential V{R) and the energy conditions EC1= p + Pr + ^Pt, EC2= p + Pr and 
EC3= /9 + pt, for 7 = -1, a; = 0.1, b = 0.000001 and rUc = 0.000866559. Case A 

can deduce from Table III is that we can have black hole formation even with an interior 
phantom energy for any given 7. The results obtained in this paper further confirm our 
previous conclusion: even though the existence of gravastars cannot he completely excluded 
in these dynamical models, our results do indicate that, even if gravastars indeed exist, they 
do not exclude the existence of black holes. 
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FIG. 3: The potential V{R) and the energy conditions EC1= p + Pr + ^Pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 0.1, b = 0.000001 and rUc = 0.0087291. Case A 




FIG. 4: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = 7/4, oj = 0.1, b = 0.000001 and rUc = 0.84374. Case B 
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FIG. 5: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 0.1, b = 0.000001 and rric = 0.51205. Case C 



y=-1.0, w=0.1 , b=0.000002 , m =0.001 2876042 



-0.5[ 



ni^=0.001 
1113=0.002 



1 1 1 1 1 




123456789 
R 



0.27 



xlO"' 



.0 



o 
> 



c 



0.2U 



w=0.1,b=0.000002 



EC1 
EC2 

EC3«EC1 



10 



FIG. 6: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = -1, a; = 0.1, b = 0.000002 and rUc = 0.0012876042. Case A 
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FIG. 7: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 0.1, b = 0.000002 and rric = 0.0115153. Case A 




FIG. 8: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 7/4, co = 0.1, b = 0.000002 and rric = 0.843415. Case B 
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FIG. 9: The potential V{R) and the energy conditions EC1= p + Pr + ^Pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 0.1, b = 0.000002 and rric = 0.51211. Case C 




FIG. 10: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = -1, a; = 0.1, b = 0.0001 and rric = 0.0120041542. Case A 
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FIG. 11: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 0.1, b = 0.0001 and rUc = 0.0546980. Case A 




FIG. 12: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = 7/4, oj = 0.1, b = 0.0001 and rUc = 0.84269. Case B 
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FIG. 13: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 0.1, b = 0.0001 and rUc = 0.51206. Case C 




FIG. 14: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = -1, a; = 0.1, b = 0.05 and rric = 0.3271379593. Case A 
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FIG. 15: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 0.1, b = 0.05 and rUc = 0.4923412. Case A 




FIG. 16: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = 7/4, to = 0.1, b = 0.05 and rric = 0.67527. Case B 
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FIG. 17: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 0.1, b = 0.05 and rUc = 0.50247. Case C 




FIG. 18: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, ioT = -1, uj = 0.2, b = 0.05 and rric = 0.3357901746. Case A 
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FIG. 19: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 0.2, b = 0.05 and rUc = 0.499554. Case A 




FIG. 20: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = 7/4, oj = 0.2, b = 0.05 and rUc = 0.66938. Case B 
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FIG. 21: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 0.2, b = 0.05 and rUc = 0.50195. Case C 




FIG. 22: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, ioT = -1, uj = 1.5, b = 0.01 and rric = 0.2097639045. Case G 
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FIG. 23: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 1.5, b = 0.01 and rUc = 0.3775. Case G 




FIG. 24: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + Pt, for 7 = 7/4, u = 1.5, b = 0.01 and rUc = 0.75153. Case H 
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FIG. 25: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 1.5, b = 0.01 and rUc = 0.50833. Case I 




FIG. 26: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= /9 + Pt, for 7 = -1, a; = 3, 6 = 0.01 and nic = 0.2552945103. Case G 
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FIG. 27: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, w = 3, 6 = 0.01 and rric = 0.426756. Case G 
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FIG. 29: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, w = 3, 6 = 0.01 and rric = 0.50661. Case I 
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FIG. 30: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= /9 + pt, for 7 = -1, a; = 5, 6 = 0.01 and rric = 0.3033168157. Case G 




FIG. 31: The potential V{R) and the energy conditions EC1= p + Pr + "^Pt-, EC2= p + Pr and 
EC3= /J + for 7 = 0, a; = 5, 6 = 0.01 and rUc = 0.472267. Case G 
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FIG. 32: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 7/4, u = 5, b = 0.01 and rric = 0.68998. Case H 




FIG. 33: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + for 7 = 3, a; = 5, 6 = 0.01 and rric = 0.50436. Case I 
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FIG. 34: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= /9 + pt, for 7 = -1, a; = 10, b = 0.000001 and rUc = 0.0025477672. Case G 
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FIG. 35: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = 10, b = 0.000001 and rric = 0.018559. Case G 
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FIG. 36: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 7/4, uj = 10, b = 0.000001 and rric = 0.84370. Case H 
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FIG. 37: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = 10, b = 0.000001 and rric = 0.51208. Case I 
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FIG. 38: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 3, a; = -1.5, b = 0.01 and rric = 0.51211. Case K 
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FIG. 39: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 7/4, lo = -1.5, b = 0.01 and nic = 0.83759. Case L 
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FIG. 40: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = -1, a; = -1.5, b = 0.01 and rric = 0.50574. Case J 
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FIG. 42: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = -1.5, b = 0.01 and rUc = 0.51735. Case J 




FIG. 43: The potential V{R) and the energy conditions EC1= p + Pr + 2pt, EC2= p + Pr and 
EC3= p + pt, for 7 = 0, a; = -1.5, b = 0.01 and rUc = 0.150217075. Case J 
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